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We derive an analytic approximation for the concurrence of weakly mixed bipartite quantum 
states - typical objects in state of the art experiments. This approximation is shown to be a lower 
bound of the concurrence of arbitrary states. 
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Entangled states constitute one of the most fundamen- 
tal differences between quantum and classical mechanics. 
The objects of prime interest in most experiments are 
pure entangled states. Therefore, experimental setups 
are designed such that unavoidable environment coupling 
is minimised, in order to keep the system state as pure 
as possible. On the theoretical side, most efforts are con- 
cerned with arbitrary mixed states - so far with limited 
success. Hardly any entanglement measure can be calcu- 
lated for systems larger than bipartite two-level systems 
- the smallest ones of interest. 

It is astonishing that theory did not proceed more 
closely along the lines indicated by experiments, although 
approximations inspired by experimental facts are fre- 
quently - and succesfuUy - used in other branches of 
physics. In this letter we present an estimate for the en- 
tanglement of almost pure states - an approach naturally 
emerging from experimental reality, where the evolution 
into a mixed state - due to non-unitary dynamics - oc- 
curs on a time scale much larger than any other relevant 
experimental time scale Q, |^ |^. Therefore, notwith- 
standing the unavoidable evolution of an initially pure 
state into a mixed one, the system state remains quasi 
pure during the period of interest, i.e. its density matrix 
has one single eigenvalue fii that is much larger than 
all the other ones. Under this condition, all terms pro- 
portional to integer powers of /i^, i > 1, are small, and 
higher order terms can be safely neglected in our sub- 
sequent treatment of concurrence, which is one of the 
generally accepted entanglement indicators. 

The concurrence c of a purejbipartite quantum state 
I*) e Til (g) 7^2 can be defined |3 as 



cW = ^Tr(|vl/)(vl/|)2 _ Tr.gl - Tr.gl + (Tr|M/)(*|)2 , 

(1) 

with the reduced density matrices gi = Tr2|5')(^'| and 
P2 = Tri|\E')(^'|; the trace over both subsystems is de- 
noted by Tr. An equivalent and widely used expression 
reads c(5') = y'z{(W\W)^^^Trg^ , where gr is either one 
of the reduced density matrices. Though, as we will see 
in the sequel, the former definition will turn out to be 
advantageous in our following generalisation for mixed 
states g^ where concurrence is defined through the con- 



vex roof 

c(f?) = inf Vkc(*,), e = Vp»|*^)(^',;|, P^>Q. 

{P.l*.)}^ i 

(2) 

Here, the infimuni is to be found among all ensembles 
{pi, l^'i)} representing g. If we start from a given de- 
composition {e.g., the eigensystem of g, g\^i) = 
i = 1, . . . , n), all ensembles can be parametrised through 
a left-unitary transformation V <E (D^^" 'sj, where pos- 
sibly N ^ n: 



N 



N 



j=i i=i 

The concurrence of a mixed state g can therefore be ex- 
pressed as k2i] 



(g) =mY,([V ^VA yt^yt]^ 



(4) 



where the tensor A is defined as 



A Irn 



Tr(|$,)(ci>i|$0$™)- 
Tn (Tr2(|$,)(<i>,|)Tr2(|<i>fe)(<i>„.|)) 
Tr2 (Tri(|$,)(ci>,|)Tri(|$fe)(<i>™|)) 
Tr(|<i>,)($,|)Tr(|<i>,)($™|)" 



(5) 



.4 is a positive, hermitian operator satisfying A^J^ — 
A%'j — A^kj . The latter symmetry is inherited from the 
symmetric definition of concurrence in eq. J^lj and is 
crucial for our further analysis. In particular, it implies 
that A can always be expressed in terms of complex sym- 
metric matrices T" as A^^ = Tf^^iTSn)* B- Whereas 
the concurrence of a pure state is characterized by a sin- 
gle matrix t :— T^, with only one non- vanishing element, 
a mixed state in general requires more than one matrix. 
The crucial idea underlying our approximation is now 
that the concurrence of a quasi pure state can still be 
well described by a single matrix r, though with more 
than just one non- vanishing element. 

This becomes clear when wc remind ourselves of the 
obvious proportionality A^J^ ^ y/JjJf^jMpm- For quasi 
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FIG. 1: Concurrence (left y-axis labels) in qpa (solid line), together with the optimized lower (dashed line) and upper (dotted 
line) bounds Q on concurrence for a 3 x 5 state g that evolves according to a random, non-unitary time evolution. Different 
degrees of mixing are indicated by the von Neumann entropy (double-dotted dashed line; right y-axis labels). The gap between 
the qpa and the upper bound increases from (a) to (c), with increasing mixing of q. Though, even for significantly mixed states 
this difference does not exceed 0.06, and the qpa thus provides a reliable estimate of concurrence in a wide parameter range. 



pure states with fii ^ fii, i > 1, this relation induces a 
natural order in terms of the small eigenvalues /i^, i > 1. 
The leading term A\l ~ /i^ (lowest order in the ^i) is 
sufficient to characterize the concurrence of the pure state 
|$i). For mixed states, however, there are also terms of 
first order, alike Ajl ^ y/JIJ, J > 1, of second order, such 
as A^l ~ y/JIjJIk, j,k > 1, and of third and fourth order. 
The approximation 



A 



jk 



with Tjk 



(6) 



is exact up to first order, and even correctly represents 
second order elements of type A^j.. Importantly, this 
quasi pure approximation (qpa) simplifies eq. Q signif- 
icantly: 



c{g)^cM^MY,\[VrV^] 



(7) 



and a closed expression for the right hand side of this 
equation is known [lo| . It can be given in terms of 
the singular values Xi of T, this is the square roots of the 
eigenvalues of the positive hermitian matrix rr^ , 



-qp 



(g) = max(Ai - ^ A^, 0) 



(8) 



i>l 



with the Xi labeled in decreasing order. 

Note that, in eq. ©, we implicitly assumed that A\l 
does not vanish, or, equivalently, that j^i) is not separa- 
ble. This does not limit, however, the range of applicabil- 
ity of our approximation. If the dominant contribution 
|$i) to a given mixed state g is separable, g will typically 
be separable anyway. 

A major advantage of our present estimate is that - 
with respect to other methods 0, im - it even further 
reduces the computational resources for the evaluation 



of the degree of entanglement of a given state: In 0, 
it was necessary to diagonalize a matrix quadratically 
larger than g, eventually followed by an optimization 
procedure. Here, we need not optimize, and only have 
to diagonalize a matrix of the size of the given statisti- 
cal operator. Consequently, significant speed-up can be 
achieved when a large number of mixed states has to be 
assessed, such as, e.g., monitoring the time evolution of 
entanglement under environment coupling [T^ . without 
significant loss of the quality of the estimation. To illus- 
trate that, let us consider a bipartite system of dimension 
3x5, coupled to an environment. The total Hamiltonian 
acting on system and bath reads H = asHgi^lb+asbHsb, 
where Hs acts on the system alone, Hst represents the 
system-bath interaction, and as and asb determine the 
strength of system dynamics and environment coupling, 
respectively. In our following example, the matrix ele- 
ments of Hg and H^b are random numbers, and the cou- 
pling constants are fixed at ag — 0.2 and asb — 0.02. 
The total state of system and bath, initially prepared in 
a pure system state tensored with the maximally mixed 
bath state, evolves under the unitary time evolution op- 
erator U = eKpiiHt), and the dynamics of the system 
state alone is obtained upon tracing over the environ- 
mental degrees of freedom, what induces decoherence, 
i.e. mixing. In fig. ^the time evolution of concurrence is 
plotted in qpa for the above parameter values, over three 
different time intervals. Also optimized lower and up- 
per 0] bounds of concurrence are shown. The degree of 
mixing of g{t) is characterized by its von Neumann en- 
tropy. While the degree of mixing is rather small in fig. 
^a), it steadily increases in figs.Hi^b) and^Jc). Nonethe- 
less, our qpa captures the actual value of the concurrence 
rather precisely over the entire time interval, even for al- 
ready significantly mixed states - i.e. beyond its initially 
anticipated range of validity. 

Finally, let us note that the quasi pure approximation 
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FIG. 2: Concurrence in qpa for the class of 3 x 3 states defined 
in eq. Q . Although Qa exhibits rather large mixing (with von 
Neumann entropy approx. equal to 1.3 . . . 1.8, for a > 0.1) 
and positive partial transpose, the qpa (solid line) detects Qa 
as entangled, and provides a rather good approximation of 
the actual value of concurrence, which is confined by upper 
(dots) and lower Q (dashed line) bounds. 



also provides a lower bound on concurrence, and, in par- 
ticular, can distinguish entangled states from separable 
ones: We have shovifn in 7] that any symmetric matrix 
E„ z.T-, with l^aP = 1 and A'j]^ = T-{TrJ*, 
defines a lower bound of concurrence. One easily verifies 
that T defined in eq. © indeed is precisely of this form, 

with Za = {T^iT/\/j2f} l^'nP- Even entangled states 



with positive partial transpose can be characterised with 
the qpa: Consider, for example, the class of states 
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Qa 



a e [0, 1] , (9) 



of a 3 X 3-system, with /3 = {l + a)/2 and 7 = \/l — 
which were introduced in 13]. Fig. |3 compares the qpa 
to upper and lower bounds Jj, ^ of concurrence, as a 
function of the parameter a. The qpa is indeed positive 
in the entire interval, i.e. the non-separability of Qa is 
detected by purely algebraic means. 

Thus, our qpa does not only provide tools for an effi- 
cient estimation of concurrence of states with moderate 
mixing, but it is even applicable to general states. We 
reckon that the concept of quasi purity which we have 
exploited here may be a remedy also for various other, so 
far virtually uncomputable entanglement measures, and 
it appears promising to check whether such approxima- 
tions have an equally large range of applicability. 

We are indebted to Andre Ricardo Ribeiro de Car- 
valho and Marek Kus for fruitful discussions, comments 
and remarks. Financial support by VolkswagenStiftung 
is gratefully acknowledged. 
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